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Abstract —The optimal performance for forward and reverse quantum Stirling cycles is studied in this 
paper. The finite time thermodynamic performance bound and optimization criteria for these cycles are 
derived. The optimal relations between the performance parameters for quantum Stirling engines, re¬ 
frigerators and heat pumps are obtained. Published by Elsevier Science Lid. 
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INTRODUCTION 

Stirling engines, refrigerators and heat pumps are thermally powered equipment with forward 
and reverse Stirling cycles. They have the advantages of high efficiency, low vibration levels, low 
noise, long service life, simple structure, etc; therefore, they have been used widely in power en¬ 
gineering, cryogenics and energy utility industries. Some authors have studied the performance 
optimization problem for the Stirling engine and refrigerator using the method of finite time 
thermodynamic analysis [1-4]. Significant results have been obtained by Ladas and Ibrahim [5], 
Wu [6], and Hu [7]. However, in the works mentioned above, the authors have all used classical 
thermodynamics to study the Stirling systems, while the quantum performance of the working 
fluid was not considered. Only recently, finite time thermodynamics has made considerable pro¬ 
gress, and the object of study now has been extended from classical thermal systems to quantum 
thermal systems [8-13]. Geva and Kosloff [9,10] proved that the optimal efficiency for a quan¬ 
tum Carnot engine at its maximum power output is the same as that derived by Curzon and 
Ahlborn [1]. Wu, Sun and Chen [14-16] obtained the optimal performance parameters for the 
quantum Carnot refrigerator and the quantum Stirling engine by solving the generalized master 
equation of an open system [17]. Based on Ref. [14] and using finite time thermodynamics and 
nonequilibrium quantum statistical theory, the performance optimizations for the quantum 
Stirling engine, refrigerator and heat pump are extended in this paper. The optimal performance 
bound and optimization criteria for the forward and reverse quantum Stirling cycles at several 
limiting cases are derived. The results provide a theoretical basis for optimal design and oper¬ 
ation of the Stirling engine, refrigerator and heat pump. Notice that, in the calculation of this 
paper, we set /?' = h/2n = 1, where It is Planck’s constant. 


QUANTUM STIRLING CYCLES 

The quantum Stirling cycle is made of two isothermal branches and connected by two regen¬ 
erating branches. The working fluid of the cycle consists of many noninteracting harmonic oscil¬ 
lators at the atomic level. The population of the oscillators (n) at thermal equilibrium can be 
obtained from the Bose-Einstein distribution [18] 
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n = l/[exp($o) - 1] (1) 

where cu>0 is the frequency of the oscillator, /l = 1/7) and T is the internal temperature of the 
working fluid. 

We assume that the cycle operates between an isothermal heat source B\ at temperature T h 
and an isothermal heat sink B 2 at temperature T c . The heat source and sink are infinitely large 
and their internal relaxation is very strong, thus the disturbance of the working fluid on the 
source or the sink can be ignored. Therefore, the source B\ and the sink B 2 are assumed to be 
in thermal equilibrium. 

The internal temperatures of the working fluid are 7) and T 2 in the two isothermal processes. 
We assume T 2 > T e (T e is the critical temperature) such that the Bose-Einstein condensation of 
the working fluid will not be considered. The second law of thermodynamics requires 
y5c > /?2 > /?i > /?h f°r the forward cycle and p 2 > fi c > Ph > Pi for the reverse cycle, where $=1/7) 
(/ = /?, c, 1, and 2). 

In the two regenerating processes, the volumes of the working fluid and the frequencies of the 
oscillator are both constants [14]. The frequencies of the oscillator in the two regenerations con¬ 
necting the isotherms are and co 2 with co l >co 2 . The forward and the reverse Stirling cycles 
are plotted in the (n, w) plane as shown in Fig. 1 and Fig. 2, respectively. 

Assuming that the internal temperature of the working fluid in the regenerating processes 
changes uniformly with respect to time [6], we have 

dT/dt = dOT'/d t) = ± l/k, A'62; 0 (2) 

where k is the time coefficient of the regenerations with constant frequencies. The signs “ + ” 
and “ — ” correspond to the processes 4-1 and 2-3, respectively. 

Heat leak, regenerator loss and the other irreversibilities, except heat resistance, are not con¬ 
sidered in this paper for simplicity. The cycle described above is the endoreversible quantum 
Stirling cycle. 


THERMODYNAMIC ANALYSIS AND CYCLE PERIOD 

The Hamiltonian of the harmonic oscillator system S, H s , is described in the following 
form [6] 


H s = coN = ma + a 


(3) 
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where N is the number operator, a + and a are the Bosonic creation and annihilation operators 
{N = a + a, [a,a + ] = 1}. 

From equation (3), the rate of change of energy of the working fluid is obtained by 

dE/dt — d(H s }/dr = (dct>/d?)n + cu(dn/dt). (4) 

Comparing equation (4) with the time derivative of the first law of thermodynamics, we obtain 
the instantaneous heat flow (Q) and the instantaneous power (P) in the following equations: 

Q =oj(dn/d?) (5) 

P — — (d<u/dt)n. (6) 

Substituting equation (1) into equation (5) and expanding the heat quantity integral to first 
order in /ijcu (i = /?, 1 or c, and 2) at the limit (n > > 1 and ftra < < 1) gives 

Q\ —m f codNl//?!, ni<n<n2 (7) 

02=*//?,. (8) 

where Q x and Q 2 are the heat exchange quantities in the isothermal processes 1-2 and 3-4, the 
constant 1 = ln(m i/co 2 ) > 0 is the natural logarithm of the frequency ratio, and m = + 1 corre¬ 
spond to the forward and reverse cycles, respectively. 

The heat exchanges of processes 2-3 and 4-1 are executed inside the engine, the refrigerator 
or the heat pump. If the regenerator is free of loss, the quantity of the regeneration is indepen¬ 
dent of the energy dissipation of the cycle. Thus, from equations (7) and (8), the work output 
from the forward cycle (or the work input required by the reverse cycle) is 

W=Q l -Q 2 = M 1/A - l//? 2 ). (9) 

In the two isothermal processes, the working fluid S becomes an open system owing to the coup¬ 
ling action between 5 and the heat source B\ (or heat sink B 2 ). Using the generalized master 
equation of an open system, and in the Heisenberg picture, one obtains the motion equation of 
the operator x in the system S [14] 

dx/d t = z(L s + r rB L S Bp B )x - ^{(ro kj ) + [x, Q k ]Qj - (co kj )“Qj[x, Q k ]} (10) 

where L s is the Liouville’s operator of the system S, T r B represents the trace of the operator of 
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the heat reservoir (B { or B 2 ), L sB represents the coupling action between the working fluid S and 
the heat reservoir (B { or B 2 ), p B is the density matrix of the heat reservoir (B\ or B 2 ), (ft>kj) + 
and (ffljk) - are constants related to the matrix elements of the density matrix p B , Qj and Q k are 
the operators of the working fluid A and [,] is the quantum Poisson’s bracket. 

Setting x = N, Q[=a + , and Q 2 = a, and neglecting the disturbance of the working fluid on 
the reservoir (L s + T rB L sB p B = 0) gives 

d(N)/dt = — (o£i — a 2 )n + a 2 (11a) 

where = [(o>i 2 ) + +(c»i 2 ) _ ] and a 2 = [(m 2 i) + +(m 2 i) _ ] are both constants. The solution of 
equation (11a) is given by 

n = n e + (n 0 - n e ) exp[-(ai - a 2 )t] (12) 

where n e is the value of n at t = 0, while 

«e = a 2 /(ai - a 2 ) (13a) 

n e is the asymptotic stationary value of n. At thermal equilibrium, « e must meet the following re¬ 
lation: 

« e = l/[exp(/ijCo) - 1]. (13b) 

Combining equation (13a) and (13b) gives oq =uexp[(l + q)pfj] and a 2 = — a exp(q^co). 
Substituting oq and a 2 into equation (11a) yields 

dn/dt = — a exp(qfSjCo){[exp(fljCo) — \]n — 1} (1 la) 

where a and q are both constants, j = h and c correspond to the isothermal processes 1-2 and 
3-4, respectively. 

Substituting equation (1) into equation (lib), expanding the time integral to the second order 
in pjco at the limit (n > >1 and < < 1) gives 

t x — m j[(dn/d(o)/(dn/dt)\d<x>mA/(P x — f} h ), c«i<(o<co 2 (14) 

t 2 = mA/(p e P 2 ) (15) 

where t\ and t 2 are the times of the isotherms 1-2 and 3-4, A [A = (l/mj—l/cn 2 )/a > 0] is a con¬ 
stant which is independent of temperature but dependent on the quantum performance of the 
working fluid, and m = ± 1 corresponds to forward and reverse cycles, respectively. 

From equation (2), the times t 3 and t 4 of the regeneration processes 2-3 and 4-1 are: 

h = U = K\/p x - l/p 2 ) (16) 

The cycle period time, t, is 

t = t\ + t 2 T t 3 + t 4 . (17a) 

Substituting equations (14) (16) into equation (17a) yields 

t = Ami 1/(7?! - fa) + l/(fa ~ fa)] + 2k(\/fa - l/fa). (17b) 


PERFORMANCE BOUND AND OPTIMIZATION CRITERIA 

Let us define the heat reservoir temperature ratio exponent r in the following form [19] 

fa/fa = (fa/faY = D’ : (18) 

where D = /?h//?c = falfa is the reservoir temperature ratio. 

The power output ( P ), the efficiency (rf) of the quantum Stirling engine, the cooling load (R), 
the COP (coefficient of performance, e) of the refrigerator, the heating load (77) and the 
COP (' F) of the heat pump can be obtained from equations (7)-(9) and (17b) and (18) in the 
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following equations: 

P —W/t — (1 D r )f(r, /? 2 )’ 

, 1 =W/Q l = 1 -D r , () < r < 1 
R =0 2 /t = D\f(r , &), 
e —Qi/W — D r /( 1 — D r ), 1< r< oo 
n=0i/t=/(r,/i 2 ), 

'P =Qx/W = 1/(1 — Z> r ), 1< r< oo 

and 

f{r, ji 2 ) = ).{Amp 2 D'[\/iP 2 D' - 4) + l/(j8 e - P 2 ) + 2*(1 - O'})-' 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 

(23) 

(24) 


(25) 


It is seen that, when r is constant, rj, e and W are all constants. The maximum value f Q (r) of f(r, 
f} 2 ) on p 2 can be obtained from equation (25). 

fo(r) - A [Am D rmr /(l + D 05 ) 2 /(D r - D) + 2k{\ - D r )] _1 . (26) 


When r is given, which means rj, e or W is given, the optimal values of P, R and /7 can be 
obtained from eqs (19), (21), (23) and (26) in the following equations, respectively: 


P 0 =(1 - D')f 0 (r), 
R 0 =D%(r), 
n 0 —fair). 


(27) 

(28) 
(29) 


P 0 =Ar][A(l + D°- 5 ) 2 ( 1 - »,)/(1 - r] - D) + 2k ,]]- 1 , (30) 

R 0 =A[A( 1 + i) 0 ' 5 ) 2 (l + c)/(D + cD - e) + 2 k/ef', (31) 

n o =A[A( 1 + D °- 5 ) 2 ('P - 1)/(1 + VD - T) + 2 k/'V]- 1 . (32) 


Equations (30)-(32) are the main optimization relations between the power output and the effi¬ 
ciency of the engine, between the cooling load and the COP of the refrigerator, and between the 
heating load and the COP of the heat pump, respectively. 

Using the extreme value theory of functions, the maximum values of power output, cooling 
load, and heating load can be obtained from equations (27)-(29) in the following equations: 


P m —A{2k + A[{\ + Z) 0 5 )(l - D 0 ' 5 ] 2 }” 1 , (33) 

R m =A{[A 0 - 5 + (A/D 05 ) + (2k/D)°.5f - 2k}~\ (34) 

ll m —A[2k(\ - D) - A( 1 + D 0 ' 5 ) 2 ]' 1 . (35) 

The extreme points corresponding to P m , R m , and 77 m are: 

mi =0.5, (36) 

r ml =(ln Dr' \n[D(2kf 5 /(A 05 + (2k) 05 + (DA)°- 5 )\, (37) 

r m j =(ln D)~ l In [D - (D + D°- 5 )(A/2k) 05 ]. (38) 

Substituting r ml , r m2 and r m3 into equations (20), (22) and (24) gives 

)/ml =l-^ = , ?CA , (39) 

e m —D(2k) 0 ' 5 /[(A) 0 ' 5 + (2kf 5 + (DA) 05 - D(2Akf 25 ], (40) 

=[!-/) + (D + D°- 5 )(A/2k)°- 5 ]-\ (41) 
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rj m , e m , and are finite time thermodynamic performance bounds for the quantum Stirling 
engine, refrigerator and heat pump, respectively. 

For the forward cycle, when r < r ml , P and ;; both decrease with the increase of r. When 
r = r ml , P = P m and r/ = i] m . When r ml <r < 1, P decreases but rj increases with increase of r. 
When r = 1, f] = ?/ c (i/ c =l — D = 1 —/t h /j8 c i s the reversible Carnot cycle efficiency) and 
P = 0. Therefore, the finite time thermodynamic optimization criterion for the quantum Stirling 
engine should be in the following ranges: 

r <r m i < 1 , ( 42 ) 

and 

>1ca < n < t] c . ( 43 ) 

Using similar approaches as those for the heat engine, we obtain the optimization criteria of the 
quantum Stirling refrigerator and heat pump in the following equations: 

1 < r < r m2 , ( 44 ) 


and 


1 < r < r m3 , 


(45) 


respectively. 


DISCUSSION 

In equations (27) (41), k represents the effect of the regenerating times (? 3 and t 4 ) on the opti¬ 
mum performance for the quantum Stirling cycles. Assuming the regenerating times in equations 
(25) and (26) are negligible (k = 0), we have 

[f(r, P 2 ))k= o =A{zfm/l 2J D r [l/03 2 Z> r - /J h ) + ~ MF', (46) 

[f(r)]k= o =UD r - D)[AmD\ 1 + D 0 ' 5 ) 2 ]~‘. (47) 

Define the factors of the power (b^, cooling load ((5 2 ) and heating load ((5 3 ) of the forward and 
reverse quantum Stirling cycles as 6i=P/P 0 , 52 = R/R 0 , and 8 3 = n/n o . Let us also define a tem¬ 
perature coefficient (6) as 0 = [(f$ 2 D r —f} h )\/[(p 2 D T — D)]. From equations (19), (21), (23) and (27)- 
(29) with k = 0, we have 

< 5 ! = < 5 , = <53 = [/(>% P 2 Wo(r)] k= o = (1 + D°- 5 ) 2 [[/6 + D/( 1 - 0 )]~ 1 . ( 48 ) 

Equation (48) indicates that the three factors are equal when k = 0. These factors have the 
maximum values <5i m = <5 2m = (5 3m = 1 when 6 = 0 O =1/(1 + D 05 ). Thus, when the regenerating 
time is neglected, the optimization criterion for the quantum Stirling cycles can be unified in the 
following form: 

e - h)\/[(i%D r - D)\ = 9 0 = 1/(1 + D 0 - 5 ), 

(0 < r < 1 for forward cycle and 1 < r < oo for reverse cycle). (49) 


CONCLUSION 

Comparing equations (30)-(32) and (39) (41) with the results in Refs [6,7], we see that the 
quantum Stirling cycle is different from the classical thermodynamic Stirling cycle. The differ¬ 
ence is because of the different character between the quantum working fluid and the classical 
thermodynamic working fluid. On one hand, the coupling action between the working fluid and 
the source B\ (or the sink B 2 ) obeys the law of quantum statistical mechanics. On the other 
hand, the character of the quantum working fluid depends on the quantum parameters 
X = ln(cn 1 /a) 2 ) and A = (l/coi—l/<w 2 )/a- For the Stirling cycles with quantum working fluids, 
such as Ne, He, D 2 etc., the quantum characters must be considered. The results in this paper 
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provide a sensitivity analysis for the optimum design and evaluation of the Stirling heat engine, 
refrigerator and heat pump. 
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